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Abstract 

We study the effects of Resonant Magnetic Perturbations (RMPs) 
on turbulence, flows and confinement in the framework of resistive 
drift-wave turbulence. This work was motivated, in parts, by experi- 
ments reported at the IAEA 2010 conference [Y. Xu et al, Nucl. Fusion 
51, 062030] which showed a decrease of long-range correlations during 
the application of RMPs. We derive and apply a zero-dimensional 
predator-prey model coupling the Drift- Wave Zonal Mode system [M. 
Leconte and P.H. Diamond, Phys. Plasmas 19, 055903] to the evolu- 
tion of mean quantities. This model has both density gradient drive 
and RMP amplitude as control parameters and predicts a novel type 
of transport bifurcation in the presence of RMPs. This model allows a 
description of the full L-H transition evolution with RMPs, including 
the mean sheared flow evolution. The key results are: i) The L-I and 
I-H power thresholds both increase with RMP amplitude \bx\, the rela- 
tive increase of the L-I threshold scales as APli oa |6j;pi/~^p7^, where 

is edge collisionality and ps is the sound gyroradius. ii) RMPs are 
predicted to decrease the hysteresis between the forward and back- 
transition, iii) Taking into account the mean density evolution, the 
density profile - sustained by the particle source - has an increased 
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turbulent diffusion compared witli tlie reference case without RMPs 
which provides one possible explanation for the density pump-out ef- 
fect. 

1 Introduction 

The H-Mode - which usually exhibits ELM bursts - is the target regime of 
operation for next step fusion devices e.g. ITER [1]. The ELMs pose a 
major threat to the plasma-facing components and should be controlled. A 
means to control ELMs relies on Resonant Magnetic Perturbations (RMPs) 
[2, 3, 5, 4]. Since there is evidence of residual turbulence in H-mode, and 
the turbulence can affect the profiles and thus stability, there is no reason 
- a priori - to neglect turbulence effects in the ELM control experiments. 
In other words, turbulence persists in H-mode. Moreover, RMPs should be 
activated before the 1st ELM, as the latter can cause large damage to the 
divertor and wall. Hence, understanding RMP effect on LH power thresh- 
old is crucial for ITER [6, 7, 8]. The RMP effect on turbulence is, to our 
knowledge, not broadly considered. We emphasize that experiments clearly 
show that RMPs have an effect on turbulence [9]. Ref. [9] showed experi- 
mentally that RMPs damp GAM Zonal Flows, as revealed by the decrease 
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of Long- Range-Correlations during RMP. We have shown in Refs. ([10, 11]), 
that drift- wave turbulence levels and associated Zonal Flows are modified 
by RMPs. We used a basic representative model, i.e. an extension of the 
Hasegawa-Wakatani system to include RMP effects, neglecting the plasma 
response. We applied a modulational stability analysis to the latter system 
and coupled it to the evolution equation for turbulence energy. The resulting 
coupled Drift Wave - Zonal Mode (DW-ZM) predator-prey model exhibits 
modified Zonal Flows. The underlying physical mechanisms of this modifi- 
cation of Zonal Flows are: i) RMPs tilt the magnetic field lines, thereby in- 
ducing a quasilinear (5j x 5h) torque, ii) The projection of the RMP-induced 
quasilinear torque yields a quasilinear stress which can compete against the 
Reynolds stress drive of Zonal Flows and hence damp them. We identified 
two asymptotic regimes, analogous to the - 'adiabatic electron' and 'fluid 
electron' - regimes of the linear Hasegawa-Wakatani model, depending on 
the strength of the RMP coupling parameter, proportional to the square of 
RMP amphtude [Tab. 1]. In the regime where electron force balance is sat- 
isfied, the extended Hasegawa-Wakatani system reduces to a single equation 
analogous to the hnear Hasegawa-Mima equation. The DW-ZM model has 
two possible states, thus a bifurcation is possible for a critical value of the 
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linear drive, i.e. input power. Below threshold, the system is in a low con- 
finement regime (L-mode like), characterized by a high turbulence energy 
and no Zonal Modes. Above threshold, the system is characterized by an 
enhanced confinement, i.e. lower turbulence energy, due to the presence of 
Zonal Modes. This is similar to the reference case without RMPs. However, 
in presence of RMPs, the threshold for the bifurcation between these two 
states depends on RMP amplitude. In the weak-RMP regime, the threshold 
increases linearly with the RMP coupling parameter. For sufficiently strong 
RMPs, the system undergoes a back-transition to L-mode. Another result is 
that the drift-wave frequency is radially modulated by the zonal density per- 
turbations. As the LH transition is thought to be triggered by Zonal Flows 
[12, 13], our model predicts that RMPs increase the LH power threshold, a 
key result consistent with experiments near-transition [3, 6, 7, 8]. Physically, 
RMPs enhance (coUisional) radial transport of electrons by couphng parallel 
transport to the tilting of magnetic fiux-surfaces. The resulting cross-field 
radial fiux can compete with the cross-field transport of polarization charge 
which is the agent of guiding-center ambipolarity breaking responsible for 
Zonal Flow formation. This simple model shows that in presence of RMPs, a 
pure Drift- Wave - Zonal Flow paradigm is no longer sufficient. With RMPs, 
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the Zonal Flows, in addition to being nonlinearly coupled to primary drift- 
waves, are also linearly coupled to secondary Zonal Density perturbations 
[11], resulting in increased damping of (secondary) Zonal Flows. Taking into 
account the mean density evolution, the density profile - sustained by the 
particle source - has an increased turbulent diffusion due to the increased 
turbulence saturation level, compared to the reference case without RMPs. 
The later result provides one possible explanation for the puzzling density 
pump- out effect. This coupling mechanism - acting here on mesoscales - is 
quite universal and also acts on profile scales. More precisely, with RMPs, 
the density profile - sustained by the particle source - is linearly coupled to 
the mean flow by the RMPs. The key, bottom hue result is that, through 
this direct coupling, the mean flow shear is decreased relative to its reference 
value without RMPs. Based on the RMP effect on Zonal Modes, we include 
mean flow shear effects, thereby extending the OD model of Ref. [14] to in- 
clude RMPs. Table 2 summarizes the group of different OD predator-prey 
models. This list is not exhaustive, we focus here on OD models directly 
related to the present study. This allows to describe RMP effects on the full 
LH transition scenario, including mean E x B flow. 
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Table 1: Qualitative impact of RMPs on Zonal Flows, in two asymptotic 
regimes. 





RMP coupling parameter 


Effect 


Weak-RMP 


|Bx 12 ^11 ^ -1 


weakly-perturbed Zonal Flows 


Strong-RMP 


12 ^11 1 


Electron force balance {Er)zM — —{n)zM 



2 Model 

We consider the extended Hasegawa-Wakatani model of Refs. [10, 11], includ- 
ing resonant magnetic perturbations (RMPs), for density n and electrostatic 
potential 0: 

d 

— n + {0,n} = V||j|| (1) 
p'^Vi0 + p^{0,Vl0} = Viiiii (2) 

where Vy denotes the parallel gradient, given - neglecting the plasma response 
- by: V|| ~ V||o + with SB^/B the radial magnetic perturbation [11]. 

The model (1,2) is closed by the generalized Ohm law, which can be 
written as: j\\ — —D\\V\\{(t) — n) with D\\ the electron parallel diffusivity 
D\\ oc l/fei- RMP effects are best described using a two-fluid picture. The 
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Table 2: Group of different OD predator-prey models. Here, 'KD' stands for 
'Kim & Diamond', LD stands for 'Leconte & Diamond', and 'MSF' stands 
for 'Mean Sheared Flow'. 



Model 


Components 


RMP 


Ref. 


Diamond et al. 


DWs, Zonal Flows 


no 


[15] 


KD 2003 


DWs, Zonal Flows, MSF, pi gradient 


no 


[14] 


LD 2012 


DWs, Zonal Flows, Zonal Density 


yes 


[11] 


this article 


DWs, Zonal Flows, Zonal Density, MSF, n gradient 


yes 





basic mechanism of RMP effect that we propose is improving over time. In 
our previous modeling [10, 11], we first thought that RMPs were generating a 
mesoscale electron current leading to an RMP-induced stress which competes 
against the Reynolds stress thus increasing Zonal Flow damping. In the 
present work, we realized that the latter RMP-induced stress corresponds 
to the projection - along the unperturbed magnetic field Bq - of quasilinear 
X 5B) torques localized around rational surfaces. Note that RMPs do also 
generate a mesoscale electron current - through the quasilinear {5v x 5B) 
electric field - refered to as 'y-independent eddy currents' in the nonlinear 
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tearing mode theory [16]. However, this mesoscale current is not the main 
actor in our theory and we neglect it in the present work (it is a 3rd-order 
nonhnear effect), the main actors are the RMP-induced quasihnear torques. 
Without RMPs, the amphtude olExB Zonal Flows is set by the competition 
between mesoscale Reynolds stress drive and neoclassical damping due to ion- 
ion friction. In presence of RMPs, mesocale quasilinear torques ((5j x 5'B) 
arise - in the poloidal direction - which couple the E x B Zonal Flows to the 
(electron) zonal diamagnetic flows, thus driving Zonal Density perturbations 
at the expense oi E x B Zonal Flows. Physically, the {Sj x S'B) torques 
are equivalent to an RMP-induced quasilinear stress which competes against 
the Reynolds stress. As only E x B Zonal Flows are directly driven by 
the (ion) Reynolds stress, the Zonal Flow damping is enhanced - over its 
neoclassical value - by RMPs. This mecanism can also be understood as the 
uniformization of the streamf unction of the electron fluid — n on mesoscale, 
i.e. mesoscale electron force balance. 

The model considered here is an extension of the Drift- Wave - Zonal 
Mode (DW-ZM) model [10], to include equilibrium density (mean density) 
dynamics, as well as equilibrium E x B flow shear. In this work, we are con- 
cerned with the evolution of turbulence energy. Zonal Flow energy and mean 
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density gradient. We assume a characteristic scale Ly ~ for the profile 
components and ~ A for the zonal components (i.e. Zonal Modes). The 
RMP effect is clear from the flux-surface average (V||j||) of the parallel gra- 
dient of the parallel current, which enters both the vorticity equation and 
the continuity equation: 

Here the helical currents j\\ are evaluated in the vicinity of their respective 
resonance surfaces. 

The flrst term on the r.h.s. of Eq. (3) is the contribution from RMP-induced 
Magnetostatic Cells (MC). The latter are due to the Modes with hehcity = 
RMP-helicity. At first glance, it seems that this term is negligeable - in the 
vicinity of the resonance surface x — Q - compared to the second term on the 
r.h.s., because k\\ oc x. 

For resistive Drift- Waves, the components of the RMP-induced MCs can be 
shown to satisfy —V\\Q5n -\- (^&a;^[(0)zM — (?T')zm] = 0, 50 = 0, which leads 
to n oc and suggests that the RMP effects cancel out, as in ideal MHD. 
However, taking into account a uniform mean E x B flow uj-e shows that 
the RMP-induced MCs have a density component ~ x/iiu}-^) close to the 
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resonance surface, and ~ l/x far from the resonance surface, thus cancehng 
most of the RMP effects, except near the resonance surface, which shows that 
RMP effects are radially localized. 

As we consider here a OD model, we model the RMP effects as uniform, 
and we take the value at the resonance surface x = 0. This is equivalent 
to neglecting the effect of MCs, as in Refs. [11] and [19]. We will extend 
this model to ID somewhere else, and properly take into account the effects 
of MCs then. An extension of the ID predator-prey model of Ref. [21] to 
include RMPs is planned. The term on the r.h.s. of Eq. (3) - proportional to 
the divergence of a flux - is similar to the RMP term coupling Zonal Modes 
[11]. This is not surprising, since both Zonal Flows and the Mean Flow are 
modified by the RMPs via quasilinear (^Sj x dh^ torques. 

After some algebra, we obtain using Eq. (3), in presence of RMPs, the 
mean density evolution as: 

+ ( Aurb + -Dresid) JY ~ -^Part = [\^/ ~ V^) 

with Aurb oc e the turbulent diffusivity - e the turbulence energy - and £)resid 
represents the surviving - short wavelength - turbulence in the pedestal and 
where -Drmp is the RMP-induced electron diffusivity [11]. We see from Eq. 
(4) that RMPs directly couple the mean density evolution to the mean electric 
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potential (0). 



2.1 Zonal Flow & turbulence energy dynamics with 
RMPs 

In order to obtain the evolution of the E x B Zonal Flow energy, we need 
to express the zonal modulations of density {n)zM in terms of that of the 
poloidal flow {vy)zM- In other words, we use a slaving approximation, justi- 
fied by the fact that DW turbulence only directly drives Zonal Flows. The 
Zonal Flows Vzf = {vy)zM and turbulence energy e thus evolve according to 
[10]: 

AT 2 f-, -DrMP V ... |2 

— = 7oA'e-7NLe - 1 - C— — ae Vzf 

dt \ Aurb + -DRMP/ 



2C^|^ZFp A ^ -DrmP V IT/ |2 2 IT/ |2 -DRMpAurb i.. 

Ps = l-<^7) TT) ael^zFl -Ps/^|V^zf| -7, —7^ I^^zf 

"I \ -t^turb + -t^RMP / -t^RMP + -t^turb 

In the turbulence energy evolution Eq. (5), RMP effects appear only in the 
last term on the r.h.s., which represents the Zonal Mode shear. 
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2.2 Zero- dimensional model 

Thus, including the mean flow shear effect on turbulence, the model takes 
the form: 

de 
it 

d\V^-^' 



aoNe - af{e)e\VzF\^ - a^eV" - a^e^ (7) 
dt 



dt 
dN 
~dt 
dT 
~dt 



= bf{e)e\VzF\'-bf{e)\VzF\' (8) 
= -coeN -ciN -C2{V + N) + r (9) 
= -doeT-diT + Q (10) 



Here, the quantities af^{e), bf^{e) and b2^{e) are effective (turbulence- 
energy dependent) coefficients, respectively given by: 



af{e) = 
bf{e) = 



1-C 



Dql^ + -Drmp 
-Drmp "' ^ 



Dqlc + D 



RMP 



ai (11) 



bi (12) 



bf{e)^,^p:^^^^f:^ (13) 

and the parameter C2 is an RMP-induced friction given by C2 = ^^me- . 

The model described by Eqs. (7-10) is an extension of the model of 
Ref. [14], to include RMPs, with one additional parameter C2 and three 
modified parameters ai^(e), b^f{e) and b^{e) in this model. The parameters 
C2, a^^{^), b^^{e) and b^{e) represent the effect of RMPs. They involve the 
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RMP particle diffusivity -Drmp [H]- The model is based upon a density- 
gradient (Vn)-driven turbulence, so we tacitly assume that increasing the 
particle flux corresponds to a power ramp-up. A more complete model with 
both Vn-driven and VTj-driven turbulence is desirable, but this is beyond 
the scope of this article. 

We are now faced with the task of closing the system of equations (7-10). 
For this purpose, we seek - as in Ref. [14] - an expression for the mean flow 
shear V oc {(f))". Additionally, we also need an expression for the mean flow 
V — {(f))'. The closure is obtained by considering radial and toroidal ion 
force balance combined with mean charge balance, including the RMP effect 
on the mean electric potential. The detailed calculation of the mean electric 
potential (0) is given in Appendix Al. 

The resulting expression is: 

9i7-> 9i7-l 9i7-> 9 i 7-> 9uiOXa\^) 

pjll + D^MP PiP + D^MP PjP' + DjiMP PjP + D^MP J^l 

(14) 

where we use the following notations: 

$ = -(0)', N^-{ny>0, T^-{T,y>0 (15) 
Q;m = - Kky > 0, 5'niom(e) = ^ X ^r—\ (16) 
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Here, g'mom(e) represents the effect of turbulent momentum diffusion, 
wliere Amom = i^o/i^i is tlie ratio of turbulent to residual momentum dif- 
fusivities. Note that direct RMP effects on the toroidal fl,ow are beyond the 
scope of this model, and thus neglected here. In the absence of momentum 
drive (Ilinom = 0), expression (14) - obtained neglecting compressional effects 
(direct RMP effects on ions) - predicts that RMPs slave the mean electric 
field to the sum of three components: the mean density gradient N, the 
mean ion temperature gradient T, and the turbulence energy e. Due to the 
strong neoclassical damping of the mean poloidal flow, ^ <^ 1, and hence 
for Ilinom = 0, the mean electric field is mainly proportional to the density 
gradient N and ion temperature gradient T, consistent with experiments 
without external momentum input. In fact, the RMPs act to modify the 
density gradient through two mecanisms: i) directly via the third term on 
the r.h.s. of Eq. (9), and ii) indirectly via the damping of Zonal Flows which 
increases turbulence energy and hence enhances turbulent particle diffusion, 
the first term on the r.h.s. of Eq. (9). 

We see from expression (14) that RMPs decrease all driving components 
of the perpendicular electric field by a factor 1 + ^^^^ . In a weak-turbulence, 
H-Mode like regime, this seems to be in accord with experiments which show 
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a) 


_HR^,p/n = 0.1 

^'RMp'^' = o■2 

.— ^lRMp/^l = 1 
^'RMP'^' = 10 
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w/t^=°-i 


b) 


— .^lRMp/^l = o■2 









0.5 



1.5 



Figure 1: a) Predicted mean electric field $ = —{4>)' v.s. mean density gradi- 
ent N — — (n)' given by expression (14), for different values of the RMP cou- 
pling parameter ^^^^ — "■'^J^' and b) Stability diagram showing relative 
variation of LH power threshold: Zonal Flow energy |VzfP v.s. - normalized 
- effective growth-rate 7^^/7c°, for different values of the - normahzed - RMP 



coupling parameter 



fJ-RMP _ DrMP 



a clear decrease of the perpendicular electric field in presence of RMPs (the 
so-called shallowing of the 'E'^ well'). Note however, that RMPs also act 
to decrease the mean density gradient through particle balance, and thus - 
even for weak turbulence - the system is clearly nonlinear and may undergo 
bifurcations between different states. The mean electric field is plotted v.s. 
mean density gradient [Fig. 1]. 
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In the weak-RMP regime <^ 1, the mean electric field simplifies to: 



1 - 



$ rsj — 

-Drmp 



1-2- 



D 



RMP 



e + 



1 - 



iV- 
-Drmp 



1 - 



D 



RMP 



T 



(17) 



whereas in the strong- RMP regime ^ 1, it simplifies to: 



plaM 



1-2 



D 



RMP 



D 



RMP 



e + 



RMP 



D 



RMP 



(18) 



Note however, that in the strong-RMP regime, direct RMP effects on ion 
dynamics - beyond the scope of this article - are no longer negligeable, and 
could play an important role. 

The mean E x B fiow is simply the opposite of the mean electric field, so 
expression (14) directly yields: 



PtP + D 



plotM 



-e — 



RMP 

pIp 



pI/jl + -Drmp PsP + -Drmp 



pIi^ + -Drmp 



-T 



(19) 



To evaluate the associated mean fiow shear, we note that in the scrape-off 
layer (SOL) - a region of the plasma not accounted for in this simple model 
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- the mean electric potential is always positive, due to sheath effects. Since 
we showed that the mean potential gradient is negative (in accord with ex- 
periments) in the region where our model is valid (just inside from the SOL), 
there is clearly a strong shear in the mean electric field. Hence, the mean 
electric field shear is clearly an increasing function of the magnitude of the 
mean electric field. Therefore, we approximate this behavior by considering 
a proportionality relation between the amplitude of the mean electric field 
shear and the amplitude of the mean electric field: 

d 



dx 



|E^| oc |E^|, i.e. V (xV 



(20) 



In the weak-RMP regime, expression (19) yields: 

nrvior 









iV + 


^ -Drmp 















n 2 



for — - — ^ 1 



(21) 



Keeping only the density drive, i.e. for /(T, e, Ilmom) = in expression (21), 
we recover in the limit -Drmp — > the oc N"^ expression obtained in Ref. 
[14]. In presence of RMPs, expression (21) predicts that for increasing RMP 
amplitudes, the sensitivity of turbulence shearing to the driving gradients 
(and Reynolds stress and momentum source) decreases in the weak-RMP 
regime. In the strong-RMP regime, the dependence on the mean density 



gradient is given by 



1 - 





N + 






1 2 




^RMP 



-Drmp 



for — - — > 1 



(22) 



Interestingly, in the strong-RMP regime, the density dependence of the - 
squared - flow shear tends to the one without RMPs, since the mean electric 
field switches from being negative and close to —{(p)' ~ {n)' + Cst for weak- 
RMPs, to being positive and close to —(0)' ~ —{n)' + Cst for strong- RMPs. 
Note that at the boundary between the two regimes, i.e. at ^^^^ — 1, the 
turbulence shearing has no dependence on the mean density gradient N, due 
to the exact cancellation of the ion diamagnetic contribution by the electron 
diamagnetic contribution, as clear from the first term on the r.h.s. of Eq. 
(19). Hence, '^^^ = 1 corresponds to a minimum of the slope in the graph 
of v.s. N, or equivalently of (0)'^ v.s. —{n)'. 



3 Numerical results 

Before presenting the numerical results, we summarize the possible states 
of the model (7-10) in a table [Tab. 3]. Near the transition from the L- 
mode like state to the ZF-dominated state (I-phase like), the mean density 
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and temperature are approximately constant - they only weakly depend on 
turbulence energy -, and the nonlinear system reduces to Eqs. (7,8) with an 
effective growth-rate Y^{r,V'^). This system has two possible states: i) a 
no-ZF state corresponding to e 7^ 0, |VzfP = and ii) a ZF-dominated state 
where the ZF evolution near transition (e is slaved to Vzf) is given - in the 
weak-RMP regime - by: 



1 d\VzF\^ 



bi dt 



as 



1 + 1 + 2C 



/i 

bi 



Vzf\- 



o2, 



(23) 



Expression (23) shows that the associated bifurcation is a transcritical bi- 
furcation. The RMPs do not modify the nature of the bifurcation, but they 
do increase the threshold of the bifurcation. A plot of Zonal Flow amplitude 
I Vzf I v.s. order parameter 7^*^^ for different values of the ratio ^^^^^ is shown 
[Fig. lb]. 

The Figures [2a,b,c] show the numerical solution of Eqs. (7-10). Figure 
[2a] shows the reference case without RMPs -Drmp = 0. The fields plotted 
are: Zonal Flow energy |VzfP, turbulence energy e, mean density gradient 
N, mean temperature gradient T, and mean electric field $ = —V. The 
onset of the L-I bifurcation (from L-mode hke state to ZF-dominated state) 
determines the onset of the I — H bifurcation (from ZF-dominated state to 
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Table 3: Possible states of the model. 



State 


Turbulence energy e 


Zonal Flow amplitude \Vzf\ 


L-mode 







I-phase, no RMPs 






I-phase, weak-RMPs 


(a* + fJ'RMp)/bi 




quiescent H-mode 









H-mode like state), and hence is potentially more important. The power 
threshold for the L-I bifurcation is given in table 3. In the case without 
RMPs, Fig. [2a] shows the three regimes: L-mode like regime between P = 
and P ~ 0.5, followed by the ZF-dominated regime - which exhibits predator- 
prey oscillations - for P ~ 0.5 — 1.5, and the quiescent H-mode like regime 
for P > 1.5. In the H-mode like regime, the mean electric field is negative, 
corresponding to an 'E'^ well', and its absolute value is the sum of the density 
and temperature gradients. Fig. [2b] shows the case with RMPs with an 
RMP parameter -Drmp = 0.1. The power threshold of the L-I bifurcation 
is clearly increased compared to the case without RMPs. This is due to 
the RMP-induced friction which modifies the power-threshold as shown in 
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[Fig. lb] . Moreover, the difference between tfie tliresfiolds of the L-I and I-H 
bifurcations is increased by RMPs - i.e. the domain of the intermediate phase 
increases. As the I-H bifurcation is due to the mean flow shear stabihzation 
of turbulence [17], the increased threshold of the I-H bifurcation is due to 
RMPs decreasing both the mean density gradient through RMP-enhanced 
turbulent particle diffusion and the sensitivity of mean flow shear to the 
mean density gradient. Finally, during the L-mode like regime and the Zonal- 
Flow dominated regime (I-phase like), the mean density gradient is strongly 
decreased by the RMPs due to the RMP-enhanced turbulent diffusion. The 
latter decrease of density gradient is a possible explanation of the so-called 
'density pump-out' effect observed in experiments. However, in this simple 
model, the temperature gradient is also decreased in the L-mode and I-phase 
like regimes, due to enhanced turbulent heat diffusion, a feature not observed 
in experiments. This is likely due to the fact that our model does not take 
into account the heat channel for electrons. As a result of RMPs, the mean 
electric field is also decreased in the H-mode like regime, corresponding to a 
shallowing of the well' consistent with experiments on DIII-D [18]. For 
higher RMP amplitude [Fig. 2c], the trend is similar, i.e. the power threshold 
for the L-I transition increases further, and the power threshold for the I- 
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H transition also increases. In [Fig. 2c], the I-H transition is not shown, 
because it occurs at very high power. Bearing in mind that in experiments 
the amount of power (or fuehng) is hmited, this imphes that for a high enough 
RMP amphtude, the LH transition will not occur, and the system will instead 
remain in the I-phase like regime, where the pedestal is limited by turbulent 
diffusion. Fig. (2d) shows the RMP effect on hysteresis. During a ramp up 
and ramp down of 'power' (or fueling) F (dash-dotted line), the turbulence 
energy without RMPs (fuU-hne) shows that the L-H transition and H-L back- 
transition occur at a different 'power' threshold F, i.e. the system exhibits 
hysteresis. Note that for this hysteresis study, we choose a different set of 
parameters than in Figs. (2a-c), because in this OD model, hysteresis only 
occurs for a restricted domain in parameter space. The dashed line in Fig. 
(2d) shows that RMPs decrease the hysteresis. The relative decrease of 
hysteresis, measured by {h — ho) /ho - where h — Fdown/rup - ranges from 
~ 15% for i^RMP = 0.01 to ~ 50% for D^mp = 0.05. Our findings - namely 
the increase of power threshold and decrease of hysteresis - have unfavorable 
implications for ITER, it suggests that ELM control experiments should seek 
the minimum RMP amplitude (within a safety margin) required to suppress 
ELMs. 
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Figure 2: Evolution of e, |VzfP, N and T as a function of 'input power' 
r = O.Olt. Parameter values are Uq = Cq = do = 1, ai = 61 = 02 = 1, 
03 = 0.1, ci = di = 0.5, /i = 1 and Ilmom = 0. The mean electric field 
amplitude E± is also plotted, a) reference case without RMPs, b) case with 
RMPs, Drmp = 0.1 and c) case with RMPs, Drmp = 0.5, and d) RMP effect 
on LH-HL hysteresis: evolution of e during a 'power' ramp up and ramp 
down with (dashed-line) and without (full-line) RMPs, for a different set of 
parameters. 
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3.1 Steady-states of the model 



To get some insight into the possible bifurcations, we apply a steady-state 
analysis. The steady-states are given by: 



aQNe-af{e)e\VzY\' -a2eV'' -a^e^ = 



hf{e)e\V^FV-hf{e)\V^^\' - 
(coe + ci + C2)A^ + C2V = r 
{do€ + di)T = Q 



(24) 
(25) 
(26) 
(27) 



together with expression (19) for the mean flow V in terms of N,T,e and 
and with af (e),6f (e) and bf{e) given by Eqs. (11-13). 



1^1 



Here, we used the ansatz (20) V (x V and we chose - for simplicity - the 
coefficient of proportionality equal to 1. 

The steady-state density and temperature gradients - near-transition - 



are approximately given - for ^ ^ - by: 



N{e) 
T{e) 



N ^ — 

Cl 

f = = Cst 
di 



^ 1 A^RMP 


-1 


. 


Cl 



g _ IX 

di 



Cst (28) 
(29) 



The detailed calculation is presented in Appendix A2. 

Hence, the turbulence energy near-transition is approximately given - for 
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^ - by: 



— - 



1^2 



ZF 



e = (30) 



where 7*^^ is an effective iinear growtli rate of tlie turbulence - modified due 
to the mean sheared fiow stabilization - given by: 



= aoiV - 02 



1 - Ml ^ ^ n ^ 



■-mom 



(31) 



We see that the mean sheared flow clearly stabilizes the turbulence, 
through the second-term on the r.h.s of Eq. (31). This effective growth- 
rate is also modified by RMP effects through the factors in Eq. (31). 
Moreover, RMPs also act by modifying the nonlinear coupling term a^^(e), 
i.e. the second term on the l.h.s. of Eq. (30). This latter modification can 
have important effects, as it depends on the turbulence energy e. We now 
restrict our steady-state analysis to the weak-RMP regime, in order to keep 
analytical calculations tractable. 

The weak-RMP regime is defined by: 

-PrMP -DrMP , /oo^ 

In the weak-RMP regime, the coefficients a^^{e), lf^{e) and h^{e) reduce 
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to: 



af(e) 
bf{e) 



bf{e) 



1-2C 
1-2C 



D 



RMP 



-Drmp 



1 + 



A'RMP 



bi 



H = Cst 



(33) 
(34) 
(35) 



where /xrmp = -DrmpPs ^- Replacing af^{e) in Eq. (30) by its expression (33), 
we obtain a quadratic equation for the turbulence energy at saturation: 



r^eff 



RMP 



-e-2C^j^ ^ 



(36) 



The physical - positive - solution is given in the weak-RMP regime by: 



7' 



eflf 



as 



1-2C 



as D 



RMP 



|^ZF|^ for ^ 



< 1 and ;^|V^zf|^ < 1 



(37) 



Apart from the solution Vzf = 0, the Zonal Flow Eq. (25) at saturation 
reduces - in the weak-RMP regime - to: 

rPs^l \ /^RMP 



1+ 1 + 2C^ 



^ , ^=0, for ^«i (38) 

Combining Eqs. (37) and (38), we obtain a logistic equation for Zonal Flow 

energy |VzfP which - in the weak-RMP regime - reduces to: 



eflf 



2_ 

ai 



1+11 + 2^^'^^^^^^^ 



DqlJ 



bi 



1-2C 



as D 



RMP 



QL 



ai 
as 



zf 







(39) 
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For the interested reader, detailed analysis is given in Appendix A2. The 
nonlinear system of Eqs. (37,39) has two possible states: i) a no-ZF state 
(L-mode) corresponding to e = 7^^/03, Vzf — and ii) a ZF-dominated state 
(I-phase) with e 7^ 0, |Vzf| 7^ 0. The effective power threshold 7^^ for the L-I 
bifurcation is given by: 



4 Discussion 

In this work we coupled the Drift- Wave Zonal Mode predator-prey model to 
the evolution of mean quantities, i.e. mean density, mean ion temperature 
and mean electric field, including direct RMP effects on mean density and 
mean electric field. The resulting predator-prey model - an extension of 
Ref. [14] - exhibits a higher power threshold than the reference case without 
RMPs, as well as a decrease of the mean density gradient, reminiscent of the 
puzzling 'density pump-out effect'. Our model also shows a shallowing of the 
Ej. well in the H-mode like regime, consistent with experiments, e.g. fig. 6b 
of Ref. [18]. A density pump-out mecanism in L-mode was presented in Ref. 
[19]. However that work neglected turbulence effects which are shown here 




(40) 
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to play an important role. Our model suggests that the density pump-out 
and the shallowing of the ^E^ well are two consequences of the same effect: 
the radial diffusion of electrons due to RMP-induced tilt of the magnetic 
field lines combined with collisions (Dy oc i^^^), else there is no irreversibility. 
Our model shows that, in the strong-RMP regime, the 'i^^' well can even 
become positive. We note however that this result should be taken with care 
as, in the strong-RMP regime, RMP effects on ion dynamics - beyond the 
scope of this article - are no longer negligeable, and could play an important 
role. We also note that although we showed that RMPs damp Zonal Flows 
and decrease the mean E x B flow shear, the mean flow itself can screen 
the RMPs [20] . This effect - beyond the scope of this work - could moderate 
the flow damping effect. There are limitations of our model. First, it is 
zero-dimensional. An extension to ID, based on the model of Ref. [21] is 
under way. Second, it neglects the plasma response, which was shown to be 
important, e.g. in Ref. [22]. Third, our model neglects direct RMP effects 
on toroidal flows. Moreover, the numerical results presented here focus - for 
the sake of simplicity - on the zero torque regime. Note however, that - in 
our model - external torque simply decreases the effective growth rate 7^^ 
as a result of increased shearing. Finally, we note that our model focuses on 



29 



zero-frequency Zonal Flows, rather than Geodesic acoustic Modes (GAM). 
A heuristic model for GAMs - including the curvature coupling can be found 
in the Appendix of Ref. [11]. 

5 Conclusion 

We investigated, in this work, RMP effects on the L-H transition. Here are 
the main results: i) The density profile - sustained by the particle source - 
has an increased turbulent diffusion compared to the reference case without 
RMPs, one possible explanation for density pump-out. ii) RMPs decrease 
the sensitivity of mean flow shear to the driving density gradient, resulting 
in a shallowing of the Er well, iii) As a result, RMPs are shown to increase 
the power threshold for both the L-I and the I-H transitions, although the 
latter effect seems stronger, and will be investigated deeper in the future. 
For reference, we summarize the scalings of the relative change in power 
threshold for the - analytically tractable - LI transition [Tab. 4]. Bearing 
in mind that in experiments the amount of power (or fueling) is limited, 
this implies that for a high enough RMP amplitude, the LH transition will 
not occur, and the system will instead remain in an I-phase like regime. 
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where the pedestal is limited, iv) RMPs can amphfy the hysteresis between 
L — H and H — L transitions. Main results on hysteresis are summarized 
[Tab. 5]. The latter findings, namely the increase of power threshold and 
decrease of hysteresis have unfavorable implications for ITER, and set a 
constraint for ELM control experiments. As RMPs should be turned on before 
the LH transition to avoid even the first large ELM, it suggests that ELM 
control experiments should not only seek an RMP amplitude compatible 
with ELM suppression, but the minimum RMP amplitude (within a safety 
margin) required to suppress ELMs. Otherwise, an unnecessary expanditure 
of power for the L-H transition will result. 

Table 4: Scalings of the relative change in the L-1 power threshold as a 
function of experimental quantities. 



RMP amplitude 


edge collisionality 


sound gyroradius 


APli<x 


B 


2 


^Pli (X 4 


^Pli (X 4 
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Table 5: RMPs decrease the LH-HL hysteresis. 



^11 2 

RMP coupling parameter — ^ 





0.02 


0.04 


0.06 


decrease in hysteresis strength (%) 





15 


40 
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Appendix Al: RMP effects on mean electric 
field 

We start from the force balance for ions: 
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and for electrons: 



(42) 



where we neglected electron inertia, and we included a (NBI-induced) volume 
force of the form: 



(43) 



Here, we depart from the usual toroidal-field approximation by consider- 



ing an axisymmetric magnetic field B not purely toroidal: 



B 



B 



(44) 



with 1^1 < 1. 



Electron Force Balance Eq. (42) yields: 



1 



B 

B 



enori 

X 



B 



-V||(/>H V||Pe 

eno 

_ 5 X 

B en^B 



(45) 
(46) 



We now solve Ion Force Balance Eq. (41) by exploiting the fact that the 
l.h.s. of Eq. (41) is small (reflected by the FLR effect ps). 
At order 0, we obtain: 



'1? 



B B V^pi 

— X H X 

B B B enoB 



(47) 
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where the subscript ± states that the flow is perpendicular to the total ax- 
isymmetric magnetic field B given by expression (44). 

Note that the two components of (46) and (50) are the usual E x B drift 
and electron/ion diamagnetic drift, respectively. However, here they have 
components due to both the toroidal magnetic field and poloidal magnetic 
field, e.g for the E x B drift.: 



VjE = = 6;, X V_l0 + —^ey X V_l0 



(48) 



Up to order 1 we obtain: 




CsB L 



X ViE ■ VVjE 



(49) 



where is given by: 




B 



Vx0 



(50) 



B 



We now apply Mean Charge Balance: 




(51) 



I.e.: 





) 







(52) 
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where the RMP-induced flux ^ bxj\\ y is given - in a quasihnear approximation 
-by: 



Ml) = 



d{(t)) d{pe 



dx dx 



(53) 



Replacing the velocities by their expression, we obtain: 



B ^ B ^ p,B 

-B X V±Pi V±Pe + eno— - x 

B B Cg B 



d 

dx 



(54) 



Now, we note that - since we neglect curvature effects - the divergence of 
both the ion and electron diamagnetic drifts vanishes, and Eq. (54) reduces 



to: 



eno— ( Vj 



B 

~B 



X 



+ 



d 

dx 



(Mi) 







(55) 



We note the following identity valid for any 2D vector field u and the 
(divergence-free) magnetic field B: 



V± • [B X u] = -B ■ (V± X u) 



(56) 



We also use the following approximation: 



V I X 



VjE • VVifi] ~ VjE • V_L V_L X Vjej 



(57) 



Using identity (56) and approximation (57), the mean charge balance (55) 
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can be written: 

-{viE ■ VxO) + = (58) 

where we defined the (normahzed) vorticity 

= -p^vi0 (59) 

Eq. (58) can be interpreted as a Mean Vorticity Equation in stationary state. 
We now use the Taylor identity: 

ViE • V±h^ = ^(viExViEy^ (60) 

Using the Taylor identity (60), we obtain the Mean Charge Balance as: 

-p'sQ^i&y) + (^-^'ll) - P'f^i'^^y) = (61) 
where we integrated radially and added a neoclassical flow damping term. 



en 2^ 



V± X ViE 



Note that the mean charge balance (61) has two important consequences, 
which are better understood in terms of the mean polarization charge (pro- 
portional to the E X B Reynolds stress) . 

i) First without RMPs, mean charge balance implies that the mean poloidal 
flow is set by the competition between the mean Reynolds stress drive and 
the neoclassical flow damping (/x). 
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ii) In presence of RMPs, the RMP-induced quasilinear flux (b^jw) - due to 
the poloidal component of the quasihnear (5j x S'B) torque - can compete 
against the mean E x B Reynolds stress, thus decreasing the mean poloidal 
flow. 

Equation (61) clearly shows that, in presence of RMPs, the neoclassical 
damping of the flow plays a fundamental role. In fact, RMPs (mainly acting 
on electrons), mediated by the neoclassical damping (acting on ions), can 
modify the (ion) mean poloidal flow, and thereby modify the mean E x B 
flow through (ion) radial force balance. 

Now, from the evolution equation for ions Eq. (41), we obtain the (ion) 
radial force balance and (ion) toroidal force balance: 

= + (62) 

^(^^^^^||) = '^residg^iViz) + (Fz) (63) 

Here we used the trivial identity: {vi\\) = {viz), and we added a diffusion 
term due to residual short- wavelength turbulence (t'resid)- 
We neglected the direct effects of RMPs on the ion dynamics, here, supposed 
to be small compared to RMPs effects on electrons, due to fast electron 
streaming along the field lines. 

Combining the (ion) radial and toroidal force balance Eqs. (62,63) with 
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the mean charge balance (61), assuming constant electron temperature — 
rprei _ Qg^^ obtain: 

-(7^.)^W-(^>^(7^.)-^('^) + W-§(^..) = (64) 



d_ 
dx 



(i^turb + i^resid) ^ {Viz) 



(F.) (65) 



where (F^) denotes the toroidal component of the volume force F, and the 
RMP-induced flux {bxj\\) is given by: 



{bxj\\) = --Drmp 



d TJ^' d 



(67) 



Here, ^^rmp = -D|||6a;P is the RMP-induced electron diffusivity [11], and we 
approximate the angular momentum flux (vfjji\\ ^ by a turbulent momentum 
diffusion. For simplicity, we neglect non-diffusive terms in the ion channel, 
since we focus on RMP effects on electrons. 

To obtain Eq. (64), we used the following approximation for the mean ion 
pressure: (pi) ~ {n){Ti). 

Our model consists of the three equations (64,65,66), together with Eq. 
(67). This model determines the mean perpendicular electric field (0)' in 
terms of mean ion temperature gradient (Ti)', mean density gradient (n)', 
turbulence energy e and torque (Fz). 
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In order to couple the present model to the Drift Wave - Zonal Mode (DW- 
ZM) model, we consider profiles with a characteristic gradient scalelength L^. 



ld_ 
e dx 



J_d_ 



(0) 



-1 



(68) 



where L„ is the usual density-gradient scalelength defined as ^ = — (1/ {n))d{n) / dx. 
Using this ansatz, the model (64,65,66) together with Eq. (67) reduces 



to: 



-N' - T' + Vy- -^V, = 

B 



^resid 



(69) 
(70) 



D 



RMP 



+ ^6 (71) 
A* 



where e = ^' = (^i)'' ^' = and = (0)'. and the function 

g{e) represents the turbulent diffusion effect: 

with the parameter A = yo/yi- We use a quasilinear approximation for the 
perpendicular Reynolds stress: 

(^'^^'?) ^ ^^^^^ = "Me > (73) 

Combining Eqs. (69,70,71), the mean electric field — $' is given by: 
-Drmp 



-N' - T' - $' 
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"M ^ _ -^J-rmom _ Q 
A* B i/j-esid 



(74) 



Solving Eq. (74) for — yields expression (14) given in the main text. 



Appendix A2: RMP effects on the L-I bifur- 
cation 

Near transition, the turbulence energy is slowly varying compared to Zonal 
Modes, and hence, its dynamics is effectively slaved to that of Zonal Modes, 
i.e. de/dt ~ 0. Replacing the mean flow V by its expression (19), Eq. (24) 
can be written: 



a^N — 02 



Mrmp 

M 

_j_ Mrmp ' ' \ Mrmp 



N 



(rp / \ nmom \ 



-af {€)\Vzf\ 



1 + 



2 a2 



I Mrmp _j_ Mrmp 



ai€^ = 
(75) 



The steady-state turbulence energy is given by: 



1 + 



2 02 



1-ML ^ X j-j 

I I MRMP ^ l\ MRMP I ^ 9raom{t) ^ 



1^ 



X 



—N — — 

Oi Oi 



i-^2iMP ^ . n \ 

Mat, /fTi /\ -^-^mom \ 

I I /^RMP ]^ I /^RMP 9mom{^) J 



Oi 



e = 



(76) 



In Eq. (76) remain two unknowns N and T. 
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The steady-state temperature gradient T is obtained using Eq. (27): 



T = T(e) = <7heat(e: 



Q 



(77) 



where 5'heat(e) = (1 + Atc)"^, with At = do/di the ratio of turbulent to 
residual heat diffusivities. 

The particle balance Eq. (26) can then be written - using (19) and (77) 



as: 



+C2 



1 _ Mrmp , /n 1 

N + rr—n^gheA^J-r + 



_j_ Mrmp ■ 



^ ;^ ^ mMP „ 1 ^ MRMP 



fl'mom(^) 



coe + ci + C2)N 

ni-nnr 



(78) 



After some algebra, we obtain the steady-state density gradient N as: 



TV = N{e) = ^part(e) 



r 

Cl 



1 + 



A*RMP 
1^ . 



C2 
Cl 



I \ Q I \ n-mom 

"1 i^l 



( 

1^ 



(79) 



where 5fpart(e) 



1 + I 1 + l+^'RMP ] 37 



(H-Ajve) ^,withAiv = 



|_ '^RMP 
1 + I 1 + l+MRMP ] 37 



Now, using expressions (77) and (79), we can write the turbulence energy 
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equation (76) as: 



9tmh[e) 



ai ai 



Mrmp 



^ _j_ Mrmp 



N{e) + 



1 + 



Mrmp 



(80) 



where 5'turb(e) is given by: 



5'turb(e) 



02 



1^ WIMP a 



Mrmp 

M 

2 I MRMP 



iV(6) + 



2 _j_ MRMP 



2"(e) - 5'mom(e, 



(81) 



Now we note that since all four quantities gpa.Tt{e), 5'heat(£)) 5'mom(e) 
and g'turb(e) are slowly- varying functions of the turbulence energy e near- 
transition, we have: 



5'part(e) ~ 5'heat(e) ~ 5'mom(e) ~ S'turb(e) ~ 1 



(82) 



Hence, the steady-state density and temperature gradients - near-transition 
- are approximately given - for ^ — > - by: 



A^(e) ~ N 



Cl 

Q 



1 + 



A*RMP 
1^ 



C2 
Cl 



Q n„ 



T{e) ~ T = ^ = Cst 



Cst (83) 
(84) 



Hence, the turbulence energy near-transition is approximately given - for 
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^ - by: 



— - 



1 KZF 



e = 



(85) 



where 7*^^ is an effective iinear growtfi rate of tlie turbulence - modified due 
to tfie mean slieared flow stabilization - given by: 



7 



-- aoN - 02 



_j_ Mrmp 



N + 



2 I Mrmp 



(86) 



We see that the mean sheared flow clearly stabilizes the turbulence, 
through the second-term on the r.h.s of Eq. (86). This effective growth- 
rate is also modified by RMP effects through the factors in Eq. (86). 
Moreover, RMPs also act by modifying the nonlinear coupling term af^{e), 
i.e. the second term on the l.h.s. of Eq. (85) . This latter modification can 
have important effects, as it depends on the turbulence energy e. We now 
restrict our analysis to the weak-RMP regime, to keep analytical calculations 
tractable. The weak-RMP regime is defined by: 



D 



RMP -DrmP 



< 1 



(87) 



Note that the condition (87) implies finite turbulence e 7^ 0. In the weak- 
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RMP regime, the coefficients ai^(e), bf^i^) and bf^{e) reduce to: 



bfie) 
bf{e) 



1-2C 
1-2C 



RMP 



DQLe 
-Drmp 



1 + 



DqLe 

/^RMP 



bi 



H — Cst 



(88) 
(89) 
(90) 



Replacing af^{e) in Eq. (85) by its expression (88), we obtain a quadratic 
equation for tlie turbulence energy at saturation: 



cff 



7 
ai 



ZF 



o^-^RMP,^. 12 n f 

e — 2C— IVzFf = 0, for 



D 



RMP 



D 



QL 



DQLe 



< 1 (91) 



The physical - positive - solution is approximately given by: 



cff 
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1 -2C 



as D 



RMP 
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QL 



Vzf|' + 2C 



RMP 



QL 



ai 



7' 



eff 



IV^zfI' 



eff 



r -DrMP ^ . , IT/ |2 ^ 7 

for — < 1 and \Vzf\ < — 



(92) 



After some algebra , the Zonal Flow evolution (8) reduces - in the weak-RMP 



regime - to: 



1 o?|^zf| 



bi dt 



el^zFp - 



1+ 1 + C 



P^bi \ Urmp 



ii 

bi 



'ZF 



(93) 



Combining Eqs. (92) and (93), we obtain a logistic differential equation for 
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1 rflV^ZFp 



61 dt 
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as D 
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eff 



IT^zfI' = 



(94) 



where 7^^ denotes the - RMP dependent - power threshold for the bifurcation, 
given by: 



7 



eff f /fRMp\ 

V / 



1+ 1 + 2C^ 



^, for ^«i (95) 



Note that, in addition to the standard quadratic nonhnearity |Vzf|^ oc 
E'^ - with E the ZF energy - of the logistic differential equation, Eq. (94) 
has an additional (RMP-induced) cubic nonhnearity oc E^. However, since 
the quadratic nonhnearity is already stabilizing (negative sign), the cubic 
term does not affect the nature of the bifurcation near threshold. Hence, 
RMPs do not modify the nature of the bifurcation, it remains a transcritical 
bifurcation, but they increase the power threshold. The nonlinear system of 
Eqs. (92,94) has two possible states: i) a no-ZF state (L-mode) corresponding 
to e = Y^/a^, VzF = and ii) a ZF-dominated state (I-phase) where - 
neglecting the cubic term - the ZF energy |VzfP is given by: 



'ZF 



1 + 2C 



2 



QL 



^cff _ ^eff 

ai 



eff 



for 



/^RMP ^ . J IT/ 12 ^ 7 

< 1 and VzF < 

(96) 
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with the power threshold 7^''^ given by expression (95). 

References 

[1] F. Wagner, Plasma Phys. Control. Fusion 49 Bl (2007). 

[2] T.E. Evans et al., Phys. Rev. Lett. 92 235003 (2004). 

[3] A. Kirk et al., Plasma Phys. Control. Fusion 53 065011 (2011). 

[4] W. Suttrop, T. Eich, J.C. Fuchs et al, Phys. Rev. Lett.106 225004 
(2011). 

[5] Y.M. Jeon et al., Phys. Rev. Lett. 109 035004 (2012). 

[6] P. Gohil, T.E. Evans, M.E. Fenstermacher et al., Nucl. Fusion 51 
103020 (2011). 

[7] S. M. Kaye, R. Maingi, D. Battagha et al., Nucl. Fusion 51 113019 
(2011). 

[8] F. Ryter et al., Nucl. Fusion 52 114014 (2012). 

[9] Y. Xu et al., Nucl. Fusion 51 063020 (2011). 

[10] M. Leconte and P.H. Diamond, Phys. Plasmas 18 082309 (2011). 

46 



[11] M. Leconte and P.H. Diamond, Phys. Plasmas 19 055903 (2012). 

[12] G.S. Xu et al., Phys. Rev. Lett. 107 125001 (2011). 

[13] P. Manz et al., Phys. Plasmas 19 072311 (2012). 

[14] E. Kim and P.H. Diamond, Phys. Rev. Lett. 90 185006 (2003). 

[15] P.H. Diamond et al.. Plasma Phys. Controlled Fusion 47, R35 (2005). 

[16] P.H. Rutherford, Phys. Fluids 16 1903 (1973). 

[17] M.A. Malkov and P.H. Diamond, Phys. Plasmas 16 012504 (2009). 

[18] K. H. Burrell et al.. Plasma Phys. Control. Fusion 47 B37-52 (2005). 

[19] Q. Yu and S. Gunter, Nucl. Fusion 49 062001 (2009). 

[20] R. Fitzpatrick and T.C. Hender, Phys. Fluids B 03 644 (1991). 

[21] K. Miki, P.H. Diamond, O. Gurcan et al., Phys. Plasmas 19 092306 
(2012). 

[22] F.L. Waelbroeck, I. Joseph, E. Nardon et al. Nucl. Fusion 52 074004 
(2012). 

47 



